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Collective flavor transformations in supernovae, caused by neutrino-neutrino interactions, are 
essentially a two-flavor phenomenon driven by the atmospheric mass difference and the small mixing 
angle 613. In the two-flavor approximation, the initial evolution depends logarithmically on 613 and 
the system remains trapped in an unstable fixed point for S13 = 0. However, any effect breaking 
exact Vh-Vt equivalence triggers the conversion. Such three-flavor perturbations include radiative 
corrections to weak interactions, small differences between the i^^ and i/t fluxes, or non-standard 
interactions. Therefore, extremely small values of 613 are in practice equivalent, the fate of the 
system depending only on the neutrino spectra and their mass ordering. 

PACS numbers: 14.60. Pq, 97.60.Bw 



I. INTRODUCTION 

Dense neutrino ensambles exhibit collective flavor 
transformations [ll-Q- For supernova (SN) neutrinos, 
these effects imprint intriguing features on the processed 
spectrum such as spectral splits f7l-[l9|. Collective oscil- 
lations originate from a characteristic instability in flavor 
space of the interacting dense neutrino system that partly 
swaps its flavor content to minimize its energy. Assuming 
the primary ly^ and fluxes exceed those of the other 
flavors, the dominant effect arises for inverted neutrino 
mass hierarchy. It is driven by the atmospheric mass 
difference and the mixing angle ^13 that is known to be 
small and could vanish entirely. 

The transformation arising from an instability implies 
that the processed spectrum is independent of the mixing 
angle as long as it is small (collective transformations 
in the presence of matter are suppressed for maximal 
mixing because one projects on the interaction direction, 
but this cos 26* effect is irrelevant if 61 < 1 [13, HI] ) . In 
a two- flavor treatment, d enters only as a trigger to the 
subsequent evolution, so in the SN context a very small 
9 shifts logarithmically the onset radius for collective 
transformations [1, @. In numerical studies, choosing 
6 as small as allowed by the machine precision barely 
impacts the processed spectrum, although for 9 = the 
system remains stuck in the unstable fixed-point solution 
defined by the initial conditions. 

Such a situation looks unphysical — it does not seem 
plausible that, at least in principle, one can distinguish 
between 9 being exactly zero and some arbitrarily small 
but non-zero value. One may speculate, for example, that 
quantum fluctuations could trigger the transformation 
even for = Q , noting that collective transformations 
actually preserve flavor lepton number. 

The purpose of our paper is to show that in real 
life we do not need to worry about such subtleties. If 
6*13 is sufficiently small, three-flavor effects [l^, [23l - [27j 
trigger the instability and the logarithmic ^13 dependence 
saturates at a small but non-zero value. 



Why are collective SN neutrino transformations an 
effective two-flavor phenomenon anyway? In the outer 
layers the temperature is too low to support thermal fi or 
T populations, obviating the possibility to distinguish be- 
tween Ufj_ and Ur by charged-current reactions. Ignoring 
radiative corrections, these flavors are exactly equivalent, 
allowing us to define new flavors and Vy such that 
effectively 6*23 = 0. If in addition 9i3 = 0, one of the 
new states, say i/y, becomes equivalent to 1/3, decoupling 
entirely from the other flavors: We are left with a two- 
flavor system consisting of i^e and Vx, governed by 9i2 
and the solar mass difference 5m^. 

This system can show collective transformations. How- 
ever, the solar mass hierarchy is normal, suppressing the 
dominant transformation effect if the primary fluxes show 
the usual excess of i/g and i>e- Moreover, the collective 
oscillation region is at larger radii because the solar mass 
difference is small, so multiple-split effects are more easily 
suppressed by adiabaticity violation [H, [11]. However, 
collective oscillations driven by the solar mass difference 
do modify the spectra in some scenarios [2l|. 

Once we allow for a small but non- vanishing ^13, 
collective ^ i^y transformations become possible that 
are driven by the atmospheric mass difference Am^ and 
occur in the usual region of large neutrino flux densities. 
Our main point is that for 6113 = these transformations 
are triggered by small perturbations of the exact u^-Vt 
equivalence because Ve and then no longer form an 
exact two-flavor system. Such perturbations include ra- 



diative corrections to the v^^ and Vr matter effect [28ll2 
or small v^f-v^ flux differences. The latter can be caused 
by the presence of muons in deeper layers of the SN 
core and by radiative corrections to the interaction rates, 
modifying the relative and v-^ opacities. Non-standard 
interactions can also break the v^-v^ symmetry [2^ . a 
possibility that we will not pursue here. 

We begin in Sec. |TT] with a brief discussion of the 
equations of motion. In Sec. IIIII we prove that for 
^13 = and for exact v^-Vt equivalence, collective 
oscillations driven by the atmospheric mass difference 
arc not possible, justifying the usual two-flavor treat- 
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ment. In Sec. IIVI we study concrete departures from 
^fi~^T equivalence in the limit = and show that 
collective transformations arc triggered by these effects. 
In schematic models we compare them with an equivalent 
^13 that would trigger collective transformations at the 
same onset radius. In Sec. |V] we consider a realistic SN 
and study the competition between a small ^13 and a 
small Vfj^-Vr flux difference. We conclude with a brief 
summary in Sec. IVII 



II. EQUATIONS OF MOTION 

A. Matrix form 

For our conceptual discussion it is sufficient to con- 
sider the simplest three-flavor system showing collective 
transformations. We take the neutrino ensemble to be 
homogeneous and isotropic, study its time evolution, and 
describe mixed neutrinos by matrices of densities qe for 
each energy mode E. We use an overbar to represent 
the corresponding quantities for antincutrinos. Diagonal 
entries are the usual occupation numbers, whereas off- 
diagonal entries encode phase information. The equa- 
tions of motion (EoM) are 



\Qe = [^e,Qe\ and xqe^IV^e, 



Qe\ 



(1) 



The Hamiltonian 3x3 matrix is made up of the vacuum, 
matter, and neutrino-neutrino terms 



Hvac 
E 



(2) 



Here H^^'^ 



UM2uV2£;, with U = R23R13R12 the 
neutrino mixing matrix and M = diag(TOi, m2, TO3) the 
mass matrix. We use the standard notation of Rij as the 
rotation matrix between the i and j mass eigenstates, 
with argument Bij . For antineutrinos, the vacuum Hamil- 
tonian picks up a relative minus sign (H™"^ = — H™'^), 
whereas all other pieces remain identical. For oscilla- 
tion studies, we may neglect terms proportional to the 
identity and may write in the mass basis 



Hvac 
E 



1 / Sriv 



(3) 



The solar mass-squared difference dm? > 0, whereas the 
atmospheric one /S.m? < for inverted mass hierarchy 
(IH) and Am^ > for normal hierarchy (NH). 

The matter term, due to neutrino interactions with the 
charged Icptons, is in the flavor basis 



= V2GFdiag{N,,N^,Nr)+OiGl) 



(4) 



where iVe is the net electron density (electrons minus 
positrons) and similarly for the other leptons. The 
second-order term is due to radiative corrections and can 
be non- negligible at high densities psj . The contribution 
associated to v-v interactions is 



H'^" = V2Gf J dE (qe ' qe) + 0(G|) . 



(5) 



Multi-angle effects are ignored in our isotropic system. 
Radiative corrections can be important in dense neutrino 
gases at the second order p9l |. 



B. New interaction basis e—x—y 

Since we are concerned with a system where the 
and Vt flavors are exactly or approximately equivalent, 
it is more useful to introduce new flavors x and y that 
simplify the mixing matrix (2^ 





|=rL| 


( 













(6) 



For 



Here "unmixes" i^^ and Vr with the angle ^23- 
6'i3 = 0, i/y is the mass eigenstate ^3. Henceforth the 
interaction basis is understood to be the e-x-y-basis. 

This basis is useful because it explicitly removes ^23 
from the formalism, if the Hamiltonian and the initial 
conditions do not distinguish and Vr- Naturally, 
the evolution of and D^. is independent of ^23 in this 
approximation. 



C. Expansion in Gell-Mann matrices 

The commutator structure of the equations of motion 
ensures that the trace of qe is conserved, so we may 
re-define them to be traceless by subtracting a term 
proportional to the identity matrix I. The traceless 
part can be expanded in Gell-Mann matrices with the 
expansion coefficients forming an 8-vector X. Thus one 
can project any matrix X as 



X = tr(X)(^+X.^ 



(7) 



where A is an 8-vector of A matrices. We normalize as 
|X| = 2/y/3, corresponding to the conventions Ai = A|, 
tr(A,) = 0, tr(A,Aj) = 26,^, and [A,, A^] = 2i/yfeAfe. 
Here fijk arc the SU(3) structure constants, where 



/123 

/l47 
/458 



1 

/246 = /257 = /345 
/678 = V3/2 



1156 



-h 



367 



1/2 



(8) 



are the non- vanishing values. 

The neutrino matrices of density can now be de- 
composed, as in Eq. (O, in terms of an 8-dimensional 
polarization vector 



0E = nE I 3 + • ^ 



(9) 



is the total neutrino density per unit energy interval. 
Analogous expressions pertain to antineutrinos. 
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The different parts of the Hamiltonian can also be 
similarly decomposed. The vacuum Hamiltonian is 



where uje 
I 



2 



l\n?l(2E). The "magnetic field" is 



(10) 



B 



^13 






s 

Sn 



V-^(1 + 3C,3)/ 



/ 



Cl3'5'i2 



'iCi2(3 + Ci3) 
2' 



■^Cl2'5'l3 



V 





— •Sl3'S'l2 



•fir 02 



(11) 



sin Qi 



with e^j = brr? I lS.vr? . Moreover, we use 
Cy = cos Qij , Sij — sin 20^ , and Cy = cos 20y . 

Ignoring a term proportional to identity, the matter 
term can be written as 



Adiag(l,0,eA) = A( ^ + L-^ 



(12) 



Here A = v^Gp^e is the effective MSW potential. 
For later reference we have included ^ 1, encoding 



radiative corrections or small 
leptonic "magnetic field" is 

ll\ I 

1 






1 

73 



flux differences. The 



s2 
'23 



V 





— s 



— 5'23 
5^(1+3^23)7 



(13) 



The v-v interaction term finally is 



= M 



3+°-2 



(14) 



where the effective neutrino-neutrino interaction energy 
is ^ = V2Gf(^ + ^)- Here N = N^^^ N^^ + N^^ is the 
overall neutrino density, and N for antineutrinos. The 
collective vector D is explicitly 



D 



AE 



He P_e - 7i_e P 



E 



N + N 
The EoM are then 

Pe = {+ujeB + XL + fiB) xPe 
= {-uje'B + XL + ^I^D) xPe 



(15) 

(16) 
(17) 



The 8-dimensional vector product is defined as (ax b)i = 
fijkdjbk- In this form, the problem resembles a set of 
polarization vectors P^; precessing under the influence of 
the combined magnetic fields B, L, and the mean field 
D due to all polarization vectors. 



III. EXACT v^,-v^ EQUIVALENCE 

In the approximation that nothing distinguishes be- 
tween the and Vt flavor, 2-3 mixing is physically 
irrelevant and we expect that oscillations reduce to a 
two-flavor problem. In fact, for Oi^ = 0, no collective 
effects driven by the atmospheric mass difference occur. 

To prove this point we study a simplified system con- 
sisting of two Bloch vectors, representing equal numbers 
of neutrinos and antineutrinos, with the single vacuum 
oscillation frequency a;. For ^13 = the magnetic field 
simplifies to 



B 



^/3 



68 -I- Elj (S'i2ei - Ci2e3) 



(18) 



where are unit vectors in the 8-dimensional flavor 
space. Assuming exact Vfi-Vr equivalence implies that 
ex = 0, and therefore 



L = e3 + i=e. 



(19) 



Likewise, if the initial v^^ and Vr densities are equal, the 
initial polarization vectors P = P are proportional to the 
same linear combination of 63 and eg . 

The static vectors B and L have components in the 1, 
3; and 8 directions, whereas the only dynamical compo- 
nent of H, the self-term D, develops an 62 component. 
The EoM of the D vector derives from the difference of 
Eqs. (HH) and ^ 

D = -e^u [^12(^3 + ^3) + Gi2(Pi + Pi)] 62 . (20) 

In other words, the vector H = wB + AL -I- /i D has only 
components in the 1, 2, 3 and 8 direction and thus can 
not mix Ve and Vy. 

The same conclusion is reached if we consider the EoM 
in matrix form. The part consisting of the e and x flavor 
and the y flavor form separate block matrices both for the 
Hamiltonian matrix and the matrices of densities. In the 
e-x-y basis and with 613 = 0, the third mass eigenstate 
z/3 is not admixed to the Ve and flavors. 



IV. BROKEN i^^-i^^ EQUIVALENCE 

Even in the absence of thermal ^ or r populations 
the exact v^-Vt equivalence is broken by several sub- 
leading effects that distinguish between these flavors. In 
this case, the 1/3 flavor does not fully decouple from the 
Ve~Vx system and collective transitions driven by the 
atmospheric mass difference are inevitably triggered. 

The flrst is provided by radiative corrections to the 
neutrino matter effect where charged leptons appear 
in the loop. Even in the absence of ordinary matter, 
similar radiative corrections arise for neutrino-neutrino 
interactions, although the detailed structure of the EoM 
becomes more complicated [l^. Since collective effects 
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require a large density of neutrinos, radiative corrections 
and thus the breaking of v^j^-Vt equivalence are unavoid- 
able. Finally we note that differences in the initial 
and Vr fluxes also provides the required instability. 



A. Radiative corrections to v-r matter effect 

The presence of matter (i.e. A 7^ 0) has a similar 
effect as decreasing the effective mixing angle, although 
in detail the dynamics is more complicated. In a frame 
rotating around L there is a fast-rotating transverse in- 
field that disturbs the system and triggers the evolu- 
tion Q. However, if matter effects distinguish and 
Vr, they can play a more important role in trigerring col- 
lective oscillations, particularly for a small mixing angle. 
The largest correction is for Vr where a background of 
ordinary matter with baryon density Nb has the same 
refractive effect on Vr and Vt that would be provide d by 
a density of real r leptons [Nf 2.6 x IQ-^ Nb) 
This subleading correction is parametrized as e\ to the 
usual matter effect A in Eq. ((T3)) . 

Off-diagonal terms in the Hamiltonian generated by 
e\ will mix Ve and Vy. In the limit ^13 — >■ 0, when 1/3 
would otherwise have decoupled, these terms play a role 
similar to 6'i3, and recouple to v^. We can estimate 
the effective 613 generated by these sub-leading matter 
effects by diagonalizing the Hamiltonian instantaneously 
in matter as 



Hvac 
E 



2E 



(21) 



where the M and U denote the mass and mixing matrix 
in matter. Using the standard parametrization, Eq. (|2ip 
can be solved for the parameters of U and M. We find 
that the matter induced 1-3 mixing is 



tan 2013 = 



2u + 2t^ujCi2 + {ex - 2)A + eaACs 



(22) 



where we ignore terms beyond the leading order in e^- 
This equation should be interpreted as providing the 
critical value of ^13 such that if ^13 < ^13, the role of v^^- 
Vr equivalence breaking is more important than 6*13 itself, 
and the two-flavor approximation is not valid anymore. 

To demonstrate this, we consider a toy model with one 
Bloch vector for neutrinos P and one for antineutrinos 
P with equal length. In the two-flavor case this would 
be the simple flavor pendulum without intrinsic angular 
momentum. A nonvanishing mixing angle triggers an 
exponential growth of the misalignment between the 
force direction and the initial orientation. The time 
(or distance) after which an 0(1) deviation from the 
initial orientation is achieved grows logarithmically with 
decreasing mixing angle. We define the radius at which 
there is a change of 1% in the t'e flavor content as 
the onset radius. In Fig. 1, we show the onset radius 
for this system, and how it depends on 6'i3. We use 
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FIG. 1: Onset radius as a function of Q\z for the toy model 
described in the text. We assume a small difference between 
and Vr refraction parameterized by ea ~ 10~^. 



A = 100 km-i, EA = 5 X 10-^ /i = 10 km-\ uj = \ 
km^^, and the mixing angles Q\2 = 0.6, 623 = 7r/4. 
Using Eq. (|22p for the chosen parameters, we expect 
that the matter induced mixing becomes important at 
6*13 « 10^^. This is in good agreement with what we 
flnd. The logarithmic increase of the onset radius stops 
below this critical mixing angle. 



B. Different primary and v-r fluxes 

Another way to break the exact v^-Vt equivalence 
is through an initial flux difference. Although this 
effect is inevitable it has not been studied in detail. 
Deep in a SN core, the temperature is large enough to 
support a thermal muon population, slightly modifying 
the primary fluxes. Moreover, the same radiative effects 
that create a refractive difference between z/^ and Vt 
also modify the scattering rates and the two flavors will 
have slightly different opacities and therefore different 
thermally driven fluxes. Obviously the discrete nature of 
particle emission and thermal fluctuations of the regions 
emitting the neutrinos would necessarily make the two 
spectra different. 

As a toy example we again assume two equal Bloch 
vectors P and P. The difference between the initial 
densities of i^^ and v-^ is parameterized as 



and the same for antineutrinos. 
effect, we have explicitly 



- , (23) 
Ignoring the matter 



D = wB X (P + P) 



(24) 



which dynamically generates components of D along es 
and 67 even though initially H has only components in 
the 1, 2, 3, and 8 directions. Therefore H develops a 
component along 65, leading to a mixing of Vp and Vy. 
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FIG. 2: Same as Fig. 1, now for equal i/f^ and Ur refraction, 
but for different fluxes as indicated by the values of eat. 



It is not straightforward to define an effective mixing 
angle in this case. The effect of the different fluxes for t/^ 
and Vr is to provide terms proportional to eAr/i to the Vx- 
Vy block in the Hamiltonian. These terms are themselves 
dynamical (time-dependent), and are communicated to 
the Vf^-Vy block by the mixing between i/g and v^- 
The effective mixing angle can be thought as being the 
initial misalignment of P from the Hamiltonian which 
is approximately proportional to e7v/(w + A). As this is 
a three-flavor effect, it must vanish when — > 0. We 
therefore expect 



'13 



U) - 



X 



(25) 



The logarithmic increase of ronsct with decreasing Oi^ 
saturates at ^13 approximately equal to the effective 
mixing ^13 induced by unequal v^-Vr fluxes. 

In Fig. 2, we plot ronsct for this system as a function of 
013 for different values of cat, illustrating this effect. We 
use the frequencies (oj, /i. A) = (1, 10, 100) km"-'^, and the 
mixing angles sin^ 9x2 = 0.314, and sin^ ^23 = 0.5. Using 



Eq. (|25p for the chosen parameters, we expect the flux- 
asymmetry induced mixing to become important at ^13 ^ 
eAf/(3 X 10^'^). This is in good agreement with what we 
find. The logarithmic increase of the onset radius stops 
below the estimated value of the mixing angle. 



V. REALISTIC SUPERNOVA 

We finally consider a more realistic SN example in a 
single-angle treatment. The neutrinos are assumed to be 
emitted isotropically from the neutrinosphere at R^, = 
10 km. We assume equal luminosities for all neutrino 
flavors, given by 



L = 1.2 X 10^^ erg/s , 



(26) 



and thermal spectra with average energies {E^^) — 10, 
[Ep^) = 15, and = 20 MeV. The electron density 



of the matter is the same as in [30| at t = 1 s after the 
bounce. For the neutrino mixing parameters we use 





= 2 X 10"^ eV^ 




= 8 X 10"^ eV^ 


sin^ 6*12 


= 0.31 , 


sin^ 6*23 


= 0.50 . 



(27) 



With these assumptions, we have calculated the onset 
radius for collective transformations as a function of 
6*13, assuming a flux difference ejq = 10~^ and ignoring 
radiative corrections to the matter effect. Our results are 
shown in Fig. 3. For 6*13 < 10"'^, the onset radius is not 
sensitive to 6*13, as expected from Eq. (|25p . 
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FIG. 3: Onset radius of collective oscillations for a realistic 
SN example, assuming ejv = 10~^. 

In a realistic SN, even when the v^-Vt equivalence is 
perfect, the onset radius depends only very weakly on ^13. 
The unavoidable breaking of this symmetry by radiative 
corrections and the the presence of charged muons in 
the deep SN core almost completely removes the 6113 
dependence in the three-flavor context. Of course, MSW 
transitions caused by the ordinary matter effect depend 
on 013 in the usual way. 



VI. CONCLUSIONS 

Collective oscillations are an instability-driven phe- 
nomenon. The system transits from its initial unstable 
configuration to a stable one, triggered by the influence 
of a disturbance. Usually one thinks of this disturbance 
as being provided by the small offset between the relevant 
fiavor and the propagation eigenstates, encoded into the 
mixing angle ^13. When this mixing angle is exactly 
vanishing, one would naively think that the oscillations 
do not take place. 

However, one should recognize that a system sitting on 
an unstable fixed point is bound to be disturbed, unless 
there arc symmetries that forbid all perturbations capa- 
ble of providing an initial disturbance. In the neutrino 
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oscillation context, this symmetry happens to be the ^- 
T symmetry - which is explicitly broken. Consequently, 
collective oscillations arc inevitable. This means that 
collective oscillations take place as usual even at 6*13 = 0, 
once triggered by subleading effects. 

Another fundamental point is that SN neutrino os- 
cillations are not sensitive to arbitrarily small values 
of the mixing angle. The fantastic sensitivity to an 
arbitrarily small mixing angle, as it appears in two-flavor 
analyses, disappears when one takes into account other 
sub-leading corrections. As a result, strategies outlined 
in Refs. [l^, [l^ may be useful for determination of the 
mass hierarchy if the relevant signals are observed, but 
not for determination of a non-zero ^13 itself. On the 
other hand, in principle we could determine the neutrino 
mass hierarchy even if ^13 were exactly zero — which 



might end up being our only hope if ^13 is beyond the 
reach of laboratory-based oscillation experiments. 
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